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Abstract. Let A (- B he cancellative abelian semigroups, and let R be an integral domain. 
We show that the semigroup ring R[B] can be decomposed, as an J?[A]-module, into a direct 
sum of ii[A]-submodules of the quotient ring of In the case of a finite extension 

of positive affine semigroup rings we obtain an algorithm computing the decomposition. 
When R[A] is a polynomial ring over a field we explain how to compute many ring-theoretic 
properties of R[B] in terms of this decomposition. In particular we obtain a fast algorithm 
to compute the Castelnuovo-Mumford regularity of homogeneous semigroup rings. As an 
application we confirm the Eisenbud-Goto conjecture in a range of new cases. Our algorithms 
are implemented in the Macaulay2 package Monomial Algebras. 



1. Introduction 

Let A <^ B he cancellative abelian semigroups, and let R be an integral domain. Denote 
by G{B) the group generated by B, and by R[B] the semigroup ring associated to B, that is, 
the free i?-moduIe with basis formed by the symbols t° for a € B, and multiplication given 
by the iti-bilinear extension of t"" ■ = t"'^^. Extending a restilt of Hoa and Stiickrad in [16j . 
we show that the semigroup ring R[B] can be decomposed, as an i?[A]-module, into a direct 
sum of i?[A]-submodules of R[G{A)] indexed by the elements of the factor group G{B)/G{A). 

By a positive affine semigroup we mean a finitely generated subsemigroup B C N™, for 
some ra. If ^ C C N™ are positive affine semigroups, K \s a. field, and the positive ra- 
tional cones C{A) C G{B) spanned by A and B are equal, then K[B] is a finitely generated 
i^[^]-module and we can make the decomposition above effective. In this case the number of 
submodules Ig in the decomposition is finite, and we can choose them to be ideals of 
We give an algorithm for computing the decomposition, implemented in our Macaulay2 |12j 
package MonomialAlgebras ^. 



By a simplicial semigroup, we mean a positive affine semigroup B such that G{B) is a sim- 
plicial cone. If B is simplicial and ^ is a subsemigroup generated by elements on the extremal 
rays of B, many ring-theoretic properties of K[B] such as being Gorenstein, Cohen-Macaulay, 
Buchsbaum, normal, or seminormal, can be characterized in terms of the decomposition, see 



Proposition 3.1 Using this we can provide functions to test those properties efficiently. 



Recall that any positive affine semigroup B has a unique minimal generating set called its 
Hilhert basis Hilb(i?). By a homogeneous semigroup we mean a positive affine semigroup that 
admits an N-grading where all the elements of IIilb(i?) have degree 1. 

One motivation for developing the decomposition algorithm was to have a more efficient 
algorithm to compute the Castelnuovo-Mumford regularity (see Section |4] for the definition) 
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of a homogeneous semigroup ring This invariant is often computed from a minimal 

graded free resolution of K[B] as a module over a polynomial ring in n variables, where n is 
the cardinality of Hilb(-B). The free resolution could have length n — 1, and if n is large (say 
n > 15) this computation becomes very slow. But in fact the Castelnuovo-Mumford regularity 
of K[B] can be computed from a minimal graded free resolution of K[B] as a module over any 
polynomial ring, so long as K[B] is finitely generated. For example, if A is the subsemigroup 
generated by elements of Hilb(-B) that lie on the extremal rays of B, and K[B] = (Bgig is 
a decomposition as graded [A] -modules, then the regularity of K[B] is the maximum of 



the regularities of the Ig as -modules (Proposition 4.1). Since the minimal graded free 
resolution of Ig has length at most the cardinality of Hilb(34) (equal to the dimension of K[B] 
in the simplicial case), and the decomposition can be obtained very efficiently, this method 
of computing the regularity is typically much faster. See Section |4] for timings. 

The Eisenbud-Goto conjecture gives a bound on the Castelnuovo-Mumford regularity [S3]. 
It is known to hold in relatively few cases. The efficiency of our algorithm allows us to test 



many new cases of the conjecture (Proposition 4.3). 

2. Decomposition 

If X C G{B) we write := {t^ | x S X}. 

Theorem 2.1. Let A C B be cancellative abelian semigroups, and let R he an integral domain. 
The R[A\-module R[B] is isomorphic to the direct sum of submodules Ig C i2[G(A)] indexed 
by elements g £ G := G{B)/G{A). 

Proof. We think of an element g G G as a subset of G{B). For 51 E G let 

■.= {heB\heg]. 

By construction, we have 

For each g ^ G, choose a representative hg ^ g C G{B). The module R ■ t^s is an 
submodule of R[B] and, as such, it is isomorphic to 

Ig:=R- {t^'^' \ber'g}C R[G{A)]. 

□ 

With notation as in the proof, we have 

This decomposition, together with the ring structure of R[A] and the group structure of G 
actually determines the ring structure of R[B]: if x G Ig^ and y £ Ig^ and xy = z as elements 
of i?[G(A)] then as elements in the decomposition of R[B] 

X-RlB]y=-p^;^Z€lg,+g,. 

Henceforward we assume that A C B C are positive afflne semigroups, and we 
work with monomial algebras over a field K. 

The set B^ = {x €z B \ x ^ B + (A \ {0})} is the unique minimal subset of B such that 
t^^ generates K[B] as a K[A]-module. We define Tg := {b G Ba \ b £ g}. Then Tg + A = T'g. 
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We can compute the decomposition of Theorem 2.1 if K[B] is a finitely generated 
module, or equivalently is a finite set. This finiteness (for positive affine semigroups 
^ C i?) is equivalent to the property C{A) = C{B), where C{X) denotes the positive ratio- 
nal cone spanned by X in Q"*. (Proof: if C{A) ^ C{B) we can choose an element x £ B on 
a ray of C{B) not in C{A), so nx G Ba for all n G N^. Thus, Ba is not finite. Conversely, 
if C{A) = C{B), then for all 6 G -B there exists G N"*" such that ni^b G A. To generate 
K[B] as a i^[^]-module, it suffices to take all possible sums of the multiples mb such that 
m < Uf, for all 6 in a (finite) generating set for the semigroup B.) Note that if Ba is finite, 
then G{B)/G{A) is also finite. 

From these observations we obtain Algorithm [T] computing the set Ba and the decomposi- 
tion oiK[B]. 

Algorithm 1 Decompose monomial algebra 

Input: A homogeneous ring homomorphism 

ip : K[yi, ...,yd]^ K[xi, ...,Xn\ 

of N™-graded polynomial rings over a field K with degy, = ej and degxj = bj such 
that ip{yi) is a monomial for all i and the gradings specify positive affine semigroups 
^ = (ei, . . . , ed) C 5 = . . . , 6„) C with C{A) = C{B). 
Output: An ideal Ig C K[A] and a shift hg G G{B) for each g e G := G{B)/G{A) with 

as Z™-graded -modules (with degt^ = b). 
1: Compute the set Ba = {b ^ B \ b ^ B + {A\ {0})}, and let {wi, . . . , Vr} be the monomials 
in K[B] corresponding to elements of Ba- For example, this can be done by computing 
the toric ideal Ib '■= ker(/3 associated to B, where 

ip:K[xi,...,Xn]^K[B], Xi^t''\ 

and then computing a monomial i^-basis f i , . . . , u r of 

K[xi, . . . ,a;„]/(/ij + V((yi, • . ■,yd)))- 

2: Partition the elements Vi by their class modulo G{A)^ forming the decomposition 

3: For each g £ G, choose a representative g £Tg. 
4: For each v £ Tg, choose c^j G Z such that 



5: Let Cgj := mm{cvj \ v G F^}. 
6: return 

h9-=9 + h ■■= K[A]{f-''^ \yerg}\g€G 



For f G Fg the element f" f is in K[A] because 
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is an expression with non- negative integer coefficients. Tlius, Ig is a monomial ideal of K[A\ 
and hg G G{B) for each 5 G G, as required. 

Example 2.2. Consider S = ((2, 0, 3), (4, 0, 1), (0, 2, 3), (1, 3, 1), (1, 2, 2)) C and the sub- 
semigroup A = ((2, 0, 3), (4, 0, 1), (0, 2, 3), (1, 3, 1)). We get the decomposition of Ba into 
equivalence classes Ba = {0, (2, 4, 4)} U {(1, 2, 2), (3, 6, 6)}. Choosing shifts hi = (-2, 0, -3) 
and /i2 = (—1, 2, —1) in G{B) we have 

K[B] ^ i^[^]{t(2A3)^ i(4,4,7)|(_^^) ^ i^[^]{t(2A3)^ i(4,4,7)|(_^^) 

= {xo,xixl){-hi) e {xo,xixl){-h2), 
where K[A] = K[xo, xi, X2, X3]/ (xfx2 — Xqx'^). 

Example 2.3. Using our implementation of Algorithm 1 in the Macaulay2 package MONO- 



mialAlgebras we compute the decomposition of Q[B] over Q[A] in case of Example 2.2 
il: loadPackage "MonomialAlgebras" ; 
i2: A = {{2, 0,3}, {4, 0,1}, {0,2, 3}, {1,3,1}}; 
i3: B = {{2, 0,3}, {4, 0,1}, {0,2, 3}, {1,3,1}, {1,2, 2}}; 
i4: S = QQ[x_0 .. x_4, Degrees=>B] ; 
i5: P = QQ[x_0 .. x_3, Degrees=>A] ; 
i6: f = map(S,P) ; 

i7: dc = decomposeMonomialAlgebra f 

o7: HashTable{ {0,0,0} => { ideal ( xq, xix| ), {-2,0,-3} } 

{5,0,0} => { ideal ( xq, xix| ), {-1,2,-1} }} 
18: ring first first values dc 

08: 



1 "o 3 



The keys of the hash table represent the elements of G. 

3. Ring-theoretic properties 

Recall that a positive affine semigroup B is simplicial if it spans a simplicial cone, or equiv- 
alently, there are linearly independent elements ei, . . . , ed £ B with G{B) = C({ei, . . . , e^}). 
Many ring-theoretic properties of semigroup algebras can be determined from the combina- 
torics of the semigroup; see \10\ \TE[ [19\ \2T\ 127] . Here we give characterizations in terms of 



the decomposition of Theorem 2.1 



Proposition 3.1. Let K be a field, B C a simplicial semigroup, and let A be the 

submonoid of B which is generated by linearly independent elements ei,...,ed of B with 
C{A) = C{B). Let Ba be as above, and K[B] = ®g^G ^g^~^9) output of Algorithm^ 

with respect to A <^ B using minimal generators of A. We have: 

(1) The depth of K[B] is the minimum of the depths of the ideals Ig. 

(2) K[B] is Cohen- Macaulay if and only if every ideal Ig is equal to K[A]. 

(3) K[B] is Gorenstein if and only if K[B] is Cohen- Macaulay and the set of shifts 
{hg}g,=G has exactly one maximal element with respect to < given by x < y if there is 
an element z £ B such that x -\- z = y. 

(4) K[B\ is Buchsbaum if and only if each ideal Ig is either equal to K[A\, or to the 
homogeneous maximal ideal of K[A\ and hg + b £ B for all b G Hilb(i?). 

(5) K[B] is normal if and only if for every element x in Ba there exist Ai, . . . , A,^ G Q 
with < Aj < 1 for all i such that x = Aje,. 

(6) K[B] is seminormal if and only if for every element x in Ba there exist Ai, . . . , A^; G Q 
with < Aj < 1 for all i such that x = Ylt=i ^i^i- 
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Proof. For every x G G{B) there are uniquely determined elements Af , . . . , G Q such that 
X = Yl'j=i ^j^j- Then by construction 

hg = '^'^.^^ min {X]\v£Tg} ej. 

Assertion (1) and (2) follow immediately; (2) was already mentioned in [27', Theorem 6.4]. 
Assertion (3) can be found in [27, Corollary 6.5]. 

(4) Let Ig be a proper ideal, equivalently, #Tg > 2. The ideal Ig is equal to the homogeneous 
maximal ideal of K[A] and hg+b G B for all b G Hilb(i?) if and only if T^ = {m+ei, . . . , m+e^} 
for some m with m + b & B for all b G Hilb(i?). Now the assertion follows from [101 Theo- 
rem 9] . 

(5) We set Da = {x £ G (B) \ x = Yfi=i hei,\i G Q and < A^ < 1 Mi}. The ring K[B] 
is normal if and only if i? = C{B) D G{B) by [18^ Proposition!]. We need to show that 
C{B) n G{B) CBii and only if Ba ^ Da- We have Ba ^ Da if and only if Da ^ Ba, since 
Ba has #G = #Da equivalence classes and by definition of Ba- Note that Da ^ G{B)nG{B) 
and Da Ci B CI Ba- The assertion follows from the fact that every element x G C{B) n G{B) 
can be written as x = x' + riiCi for some x' G Da and Ui G N. 

(6) We set Da ■-= {x e B \ x = Eli^iei,Ai G Q and < A^ < 1 Vf}. By [ig Propo- 
sition 5.32] and [211 Theorem 4.1.1] -fC[-B] is seminormal if and only if Ba ^ -D^, provided 
that ei,...,ed G Hilb(i?). Otherwise there is a /c G {l,...,d} with efc = e'^ + e';! and 
e'fc, el G B \ {0}. We set A' = (ei, . . . , e'^, . . . , e^) and A" = (ei, . . . , e'^', . . . , e^). _ Clearly 
C{A) = C_{A') = C{A"). We need to show that Ba Q Da if and only if Ba' Q Da'- Let 
X G Ba \ Da- If x — ^ B, then x G \ .Da'- If 2; — e'^ G 5, then x — e'^ G -Ba" \ -D^". Let 
X G -Ba' \ -D^', say x = "Ylij^kk ^j^j + ^k^'k ^^"^ -^i > some j. If j / /c, then x G i?A \ Da- 
Let J = A;; consider the element y = x + e'l — Ylj:^k''^j^j ^ ^ some rij G N such that 
^jjLk ''^j niaximal. It follows that y G Ba \ Da and we are done. □ 

Note that normality of positive affine semigroup rings can also be tested using the imple- 
mentation of normalization in the program NORMALIZ [6j. We remark that from Proposi- 



tion 3.1 it follows that every simplicial affine semigroup ring K[B] which is seminormal and 
Buchsbaum is also Cohen-Macaulay. This holds more generally for arbitrary positive affine 
semigroups by [7, Proposition 4.15]. 



Example 3.2 (Smooth Rational Monomial Curves in P ). Consider the simplicial semigroup 
B = ((a,0), (a-l,l),(l,a- l),(0,a)) C and set A = ((a, 0), (0, a)), say K[A] = K[x,y]. 
Note that we have a equivalence classes. We get 

K[B] - K[x, yf e y) (x"-^ y^) e . . . (x, y"-^), 

as -fr[x, y]-modules, where the shifts are omitted. In the decomposition each ideal of the 
form {x^,y^), 1 < i,j < a — 3 with i + j = a — 2 appears exactly once. Hence K[B] is not 
Buchsbaum for a > 4, since (x°~^,y) is a direct summand. In case a = 4 there is only one 
proper ideal I4 = {x,y) and /14 = (2,2); in fact (2,2) + Hilb(i3) C B and therefore K[B] is 
Buchsbaum. It follows immediately that K[B] is Cohen-Macaulay for a < 3, Gorenstein for 
a < 2, seminormal for a < 3, and normal for q < 3. Note that we could also decompose 
K[B] over the subring i^[^], where A = {{2a, 0), (0, 2a)) = K[x' , y'], for a = 4 we would get 

K[B]^K[x',yf'®{x',y') 

and the corresponding shift of {x',y') is again (2,2). 
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Example 3.3. Let B = ((1, 0, 0), (0, 1, 0), (0, 0, 2), (1, 0, 1), (0, 1, 1)) C N^ moreover, let 
A= ((1,0,0), (0,1,0), (0,0, 2)), say K[A] = K[x,y,z]. This example was given in [211 Ex- 
ample 6.0.2] to study the relation between seminormality and the Buchsbaum property. We 
have 

K[B]^K[A](B{x,y){-{0,0,l)), 

as Z^-graded i^[j4]-modules. Hence K[B] is not Buchsbaum, since {x, y) is not maximal; 
moreover, K[B] is seminormal, but not normal. 

Example 3.4. Consider the semigroup B = ((1, 0, 0), (0, 2, 0), (0, 0, 2), (1, 0, 1), (0, 1, 1)) C N^, 
and set A = ((1,0,0), (0, 2, 0), (0, 0, 2)). We get 

K[B] ^ K[A] e K[A]{-{1, 0, 1)) e K[A]{-{Q, 1, 1)) © K[A]{-{1, 1, 2)). 

Hence K[B] is Gorenstein, since (1,0, 1) + (0, 1, 1) = (1, 1, 2). Moreover, K[B] is not normal, 
smce (1,0,1) = (1,0,0) + 5(0,0,2), but seminormal. 

Example 3.5. We illustrate our implementation of the characterizations given in Proposi- 
tion [3]T] at Example [3]4j 

il: B = {{1,0,0}, {0,2,0}, {0,0, 2}, {1,0,1}, {0,1,1}}; 

12: isGorensteinMA B 

o2 : true 

13: isNormalMA B 

o3: false 

14: isSeminormalMA B 
o4: true 

Note that there are also commands IsCohenMacaulayMA and isBuchsbaumMA available testing 
the Cohen-Macaulay and the Buchsbaum property, respectively. 

4. Regularity 

Let be a field and lei R = K[xi, . . . ^ x„] be a standard graded polynomial ring, that is, 
degXj = 1 for alH = 1, . . . , n. Let i?+ be the homogeneous maximal ideal of i?, and let M be 
a finitely generated graded i?-module. We define the Castelnuovo-Mumford regularity regM 
of M by 

regM := max |a(i7|j^ (M)) + i | i > o| , 

where a{Hj^^{M)) := maxjn | [H'^^{M)]n / o| and a(0) = -00; H}^^{M) denotes the i-th 

local cohomology module of M with respect to Note that regM can also be computed in 
terms of the shifts in a minimal graded free resolution of M. An important application of the 
regularity is that it bounds the degrees in certain minimal Grobner bases by p!]. Thus, it is of 
interest to compute or bound the regularity of a homogeneous ideal. The following conjecture 
(Eisenbud-Goto) was made in [9]: If is algebraically closed and / is a homogeneous prime 
ideal of R then for S = R/I 

reg S < deg S — codim S. 

Here degS" denotes the degree of S and codim := dim/^ 5i — dimS* the codimension. The 
conjecture has been proved for dimension 2 by Gruson, Lazarsfeld, and Peskine [13j; for the 
Buchsbaum case by Stiickrad and Vogel [29] (see also [30]); for deg 5 < codim S" + 2 by Hoa, 
Stiickrad, and Vogel pTT]; and in characteristic zero for smooth surfaces and certain smooth 
threefolds by Lazarsfeld ^20j and Ran [25j. There is also a stronger version in which S is 
only required to be reduced and connected in codimension 1; this version has been proved in 



DECOMPOSITION OF SEMIGROUP ALGEBRAS 



7 



dimension 2 by Giaimo in [llj. For homogeneous semigroup rings of codimension 2 the con- 
jecture was proved by Peeva and Sturmfels [M]- Even in the simpHcial setting the conjecture 
is largely open, though it was proved for the isolated singularity case by Herzog and Hibi |15j : 
for the seminormal case by [22] ; and for a few other cases by \16\ [23] . 

We now focus on computing the regularity of a homogeneous semigroup ring -fC[-B]. Note 
that a positive affine semigroup B is homogeneous if and only if there is a group homo- 
morphism deg : G{B) — )■ Z with deg6 = 1 for all b G Hilb(-B). We always consider the 
i?- module structure on K[B] given by the homogeneous surjective algebra homomorphism 
R K[B],Xi I—)- where Hilb(i3) = Generalizing the results from [16], the 

regularity can be computed in terms of the decomposition of Theorem [2. 1| as follows: 



Proposition 4.1. Let K he an arbitrary field and let B C N"* be a homogeneous semigroup. 
Fix a group homomorphism deg : G{B) — )• Z with deg 6 = 1 for all b G Hilb(i?). Moreover, let 
A be a submonoid of B with Hilb(^) = {ei, . . . ,ed}, degCj = 1 for all i, and C{A) = C{B). 
Let K[B] = ®g^G ^9i~^9) output of Algorithm^ with respect to AQ B. Then 

(1) iegK[B] = max{reg/g + deghg \ g G G}; where regig denotes the regularity of the 
ideal Ig C K[A] with respect to the canonical K[xi, . . . , Xd\-module structure. 

(2) degK[B] = #G- deg K[A]. 

Proof (1) Consider the T = K[xi, . . . , x^j-module structure on K[B] which is given by T ^ 
K[A] C K[B], Xi ^ Since C{A) = C{B) we get by [5, Theorem 13.1.6] 

as Z-graded T-modules (where is the homogeneous maximal ideal of By the 

same theorem we obtain Wj^y^-^^{K[B]) = Then the assertion follows from 

K[B] ^ 0^gc^/g(-deg/ig) as Z-graded T-modules. 

(2) Follows from deglj, = deg/C[A] for ah 5 E G. □ 

Using Proposition 4.1 we obtain Algorithm [2J 



Algorithm 2 The regularity algorithm 

Input: The Hilbert basis Hilb(S) of a homogeneous semigroup B C N*" and a field K. 
Output: The Castelnuovo-Mumford regularity reg 
1; Choose a minimal subset {ei, . . . , e^} of Hilb(S) with C({ei, . . . , e^}) = C{B), and set 
A = (ei, . . .,ed). 

2: Compute the decomposition K[B] = ©^g^ /ig) over K[A] by Algorithm [ij 

3: Compute a hyperplane H = {{ti, . . . ,tm) G M"^ | X^j^i ^i^i — '^1 with c 7^ such that 

Hilb(5) C H. Define deg : M'" ^ M by deg(ti, ...,tm) = (EJli ajtj)/c. 
4: return reg K[B] = max{reg/g -|- deg kg \ g G G}. 



By Algorithm [2] the computation of reg K[B] reduces to computing minimal graded free 
resolutions of the monomial ideals Ig in K[A\ as /'^[xi, . . . , Xd]-modules. 

Example 4.2. We apply Algorithm [2] using the decomposition computed in Example 2.3 A 
resolution of / = (xq, xix^) as a T = Q[xo, xi, X2, xsj-module is 

^ T (-4) e T (-5) A T (-1) e T (-3) ^1^0 
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with 

/ 2 2 2 \ 

Y -Xo -X1X2 J ' 

hence reg/ = 4. The group homomorphism is given by deg6 = (bi + 62 + &3)/5 and therefore 
reg Q[-B] = max {4 - 1, 4 - 0} = 4. 

With respect to timings, we first focus on dimension 3 comparing our implementation of 
Algorithm [2] in the Macaulay2 package MonomialAlgebras (marked in the tables by 
MA) with other methods. Here we consider the computation of the regularity via a minimal 
graded free resolution both in Macaulay2 (M2) and Singular [S] (S). Furthermore, we 
compare with the algorithm of Bermejo and Gimenez [2]. This method does not require the 
computation of a free resolution, and is implemented in the Singular package mregular.lib 
[3] (BG-S) and the Macaulay2 package Regularity [26] (BG-M2). For comparability 
we obtain the toric ideal Ib always through the program 4ti2 [14J, which can be called 
optionally in our implementation (using [28j). We give the average computation times over 
n examples generated by the function randomSemigroup(a,d,c,nuin=>n,setSeed=>true). 
Starting with the standard random seed, this function generates n random semigroups B C N"^ 
such that 

• dim K[B] = d. 

• codim K[B] = c; that is the number of generators of -B is d + c. 

• Each generator of B has coordinate sum equal to a. 

All timings are in seconds on a single 2.7 GHz core and 4 GB of RAM. In the cases 
marked by a star at least one of the computations ran out of memory or did not finish within 
1200 seconds. Note that the computation of reg/g in step 4 of Algorithm [2] could easily be 
parallelized. This is not available in our Macaulay2 implementation so far. 

The next table shows the comparison for = Q, d = 3, a = 5, and n = 15 examples. 



c 


1 


2 


3 


4 


5 


6 


7 


8 


9 


MA 


.073 


.089 


.095 


.10 


.13 


.14 


.14 


.19 


.16 


M2 


.0084 


.0089 


.011 


.017 


.043 


.10 


.45 


2.8 


21 


S 


.0099 


.0089 


.011 


.013 


.020 


.046 


.18 


1.1 


6.8 


BG-S 


.016 


.030 


.19 


1.2 


15 


24 


59 


44 


77 


BG-M2 


.036 


.053 


.47 


1.8 


9.0 


19 


34 


39 


43 


c 


10 


11 


12 


13 


14 


15 


16 


17 


18 


MA 


.21 


.26 


.22 


.26 


.29 


.30 


.31 


.36 


.47 


M2 


180 


* 


* 


* 


* 


* 


* 


* 


* 


S 


30 


* 


* 


* 


* 


* 


* 


* 


* 


BG-S 


170 


520 


* 


* 


* 


* 


360 


460 


350 


BG-M2 


85 


150 


140 


250 


310 


290 


300 


410 


320 



For small codimension c the decomposition approach has slightly higher overhead than 
the traditional algorithms. For larger codimensions, however, both the resolution approach in 
Magaulay2 and Singular and the Bermejo-Gimenez implementation in Singular fail. The 
average computation times of the Regularity package increase significantly, whereas those 
for Algorithm [2] stay under one second. The traditional approaches become more competitive 
when considering the same setup over the finite field K = Z/101, but are still much slower 
than Algorithm [2] 
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c 


1 


2 


3 


4 


5 


6 


7 


8 


9 


MA 


.072 


.088 


.093 


.10 


.12 


.13 


.13 


.19 


.16 


M2 


.0075 


.0095 


.010 


.013 


.020 


.032 


.090 


.40 


2.8 


S 


.0067 


.010 


.011 


.015 


.023 


.041 


.16 


.99 


6.3 


BG-S 


.017 


.020 


.031 


.052 


.094 


.12 


.18 


.34 


.42 


BG-M2 


.030 


.037 


.064 


.14 


.34 


.48 


.80 


1.5 


2.0 


c 


10 


11 


12 


13 


14 


15 


16 


17 


18 


MA 


.21 


.25 


.22 


.25 


.29 


.29 


.31 


.35 


.39 


M2 


26 


* 


* 


* 


* 


* 


* 


* 


* 


S 


28 


250 


* 


* 


* 


* 


* 


* 


* 


BG-S 


.57 


.88 


.88 


1.1 


1.4 


1.5 


1.7 


2.5 


2.4 


BG-M2 


3.3 


4.4 


4.4 


6.4 


7.9 


7.8 


9.2 


12 


13 



Note that over a finite field there may not exist a homogeneous hnear transformation such 
that the initial ideal is of nested type, see for example [21 Remark 4.9]. This case is not covered 
and hence does not terminate in the implementation of the Bermejo-Gimenez algorithm in 
the Regularity package. In the standard configuration the package mregular.lib can 
handle this case, but then does not perform well over a finite field in our setup. Hence we 
use its alternative option, which takes the same approach as the Regularity package and 
applies a random homogeneous linear transformation. 

Increasing the dimension to d = 4 we compare our implementation with the most com- 
petitive one, that is, mregular.lib {K = Z/101, a = 5, n = 1). Here also the Singular 
implementation of the Bermejo-Gimenez algorithm fails: 



c 


4 


8 


12 


16 


20 


24 


28 


32 


36 


40 


44 


48 


52 


MA 


.13 


.31 


3.8 


13 


.69 


2.2 


1.7 


1.9 


1.5 


4.4 


6.0 


8.9 


13 


BG-S 


.61 


2.2 


46 


150 


380 


840 


940 


* 


* 


* 


* 


* 


* 



To illustrate the performance of Algorithm [2] we present the computation times {K 
Z/101, n = 1) of our implementation for d = 3 and various a and c: 



a\c 


4 


8 


12 


16 


20 


24 


28 


32 


36 


40 


44 


48 


52 


3 


.083 


























4 


.073 


.10 


.24 






















5 


.11 


.13 


.15 


.22 




















6 


.11 


.31 


.21 


.22 


.27 


.75 
















7 


.10 


.16 


.18 


.24 


.29 


.86 


1.0 


1.4 












8 


.11 


.22 


.26 


.31 


.35 


.54 


.67 


.85 


1.2 


3.6 








9 


.13 


.25 


.31 


.38 


.56 


.64 


.77 


.98 


1.4 


3.8 


5.7 


8.6 


13 



The following table is based on a similar setup for d = 4: 



a\c 


8 


16 


24 


32 


40 


48 


56 


64 


72 


80 


3 


.18 


.51 


















4 


.26 


.32 


.54 
















5 


.31 


13 


2.2 


1.9 


4.4 


8.9 










6 


9.6 


120 


* 


* 


3.4 


7.8 


15 


36 


66 


120 
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Obtaining the regularity via Algorithm [2] involves two main computations - decomposing 
K[B] into a direct sum of monomial ideals Ig Q K[A] via Algorithm [l] and computing a 
minimal graded free resolution for each Ig. The computation time for the first task is in- 
creasing with the codimension. On the other hand the complexity of the second task grows 
with the cardinality of Hilb(^), which tends to be small for big codimension. This explains 
the good performance of the algorithm for large codimension observed in the table above. In 
particular the simplicial case shows an impressive performance as illustrated by the following 
table for simplicial semigroups with d = 5 and a = 5 (same setup as before). The examples 
are generated by the function randomSemigroup using the option simplicial=>true. 



c 


10 


20 


30 


40 


50 


60 


70 


80 


90 


100 


110 


120 


MA 


13 


13 


17 


32 


69 


86 


110 


170 


250 


400 


650 


1000 



In case of a homogeneous semigroup ring of dimension 2 the ideals Ig are monomial ideals 
in two variables. Hence we can read off reg/^ by ordering the monomials with respect to 
the lexicographic order (see, for example, [231 Proposition 4.1]). This further improves the 
performance of the algorithm. 

Due to the good performance of Algorithm [2] we can actually do the regularity computation 
for all possible semigroups B in N*^ such that the generators have coordinate sum a for some 
a and d. This confirms the Eisenbud-Goto conjecture for some cases. 

Proposition 4.3. The regularity ofQ[B] is bounded by degQ[B] — codim Q [i?] , provided that 
the minimal generators of B in N'^ have fixed coordinate sum a for d = 3 and a < 5, for 
d = A and a < 3, as well as for d = 5 and a = 2. 

Proof. The list of all generating sets Hilb(i?) together with reg Q[-B] , deg Q[i?] , and codim 
can be found under the link given in [4J. □ 

Figure [T] depicts the values of degQ[i?] — codimQ[i?] plotted against regQ[-B] for all semi- 
groups with a = 3 and d = 4. For the same setup Figure [2] shows regQ[i?] on top of 
codimQ[i?] plotted against degQ[-B]. The line corresponds to the projection of the plane 
regQ[i?] — degQ[i?] + codimQ[i?] = 0. Figures for the remaining cases can be found at [4]. 
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Figure 1. degQ[B] - codimQ[i?] against regQ[i?] for a = 3 and d = A. 
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Figure 2. regQ[B] + codim Q[S] against degQ[S] for a = 3 and d = A. 
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